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ABSTRACT 

We investigate the possibihty that the difference between the measurements of as{Mz) 
from the hadronic branching ratio of the and the world average of other measurements 
is due to the decay of the Z'^ into quark, anti-squark, and gluino. Consequences for 
supersymmetry breaking models are discussed. 
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Recently, some of us have pointed out several features of quarkonium decay [|T| suggest- 
ing the existence of gluinos in the sub-GeV mass region. Gluinos below 1 GeV are natural 
in the mi/2 = model of soft SUSY breaking, which corresponds to the low-energy man- 
ifestation of spontaneously broken = 1 supergravity with minimal gauge kinetic term 
. In this model, the superpartners of the massless gauge bosons are themselves massless 
at tree level receiving small calculable masses in perturbation theory, and SUSY breaking 
is dominantly driven by a universal scalar squared mass mg representing the difference 
between the average squared masses of the fermions and their superpartners. Assuming 
degenerate super-heavy particles at the GUT scale, minimal SUSY unification with light 
gluinos and the world average values of the standard model parameters require at the 
one standard deviation level that 

75 GeV < mo < 270 GeV. (1) 

Non-degenerate super-heavies could loosen this prediction Assuming negligible left- 
right mixing for the light quark partners, the lower limit in Eq. (|1|) insures that none 
of the sfermions except perhaps the lightest stop quark is below half the mass. The 
Mz/^ experimental lower limit on the stop quark constrains the other parameters of the 
^1/2 = model. This entire range of mo leads naturally in the mi/2 = model to gluino 
masses in the sub- GeV region. Values of mo near the lower end of this allowed range are 
characterized by squark masses between Mz / 2 and Mz and it is the purpose of this paper 
to explore the consequences of this possibility. Contrary results from hadron colliders 
suggesting 100 GeV level lower bounds for the masses of squarks and gluinos are dependent 
on the assumption that these particles, once produced, will lead to significant amounts 
of missing transverse energy due to their ultimate decay into photinos or other weakly 
interacting stable neutralinos. The hadronization Monte Carlos leading to the quoted 
bounds are reliable at some level for gluino masses in the multi-GeV or tens of GeV range 
but are presumably not intended to be used for light gluinos. The following plausible 
scenario assuming squark masses above half the mass and a gluino mass below one 
GeV would invalidate the hadron collider limits. A squark produced in hadronic collisions 
will immediately decay into the corresponding quark plus a gluino. The decay gluino or a 
gluino independently produced in the hadronic collision would hadronize into a heavy jet 
of hadrons one of which would ultimately be the lightest hadron with a gluino constituent. 
Such a hadron would presumably be in the 1.5 to 2.5 GeV mass region as expected for gluon 
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bound states. The gluino-containing hadron would decay predominantly into a multi-pion 
final state plus a (predominantly soft) photino. Such a photino would very rarely carry off 
a significant amount of transverse energy. Thus if the gluinos are in the sub-GeV region 
the major effect in hadronic collisions will be only a somewhat higher apparent value of 
the strong coupling constant with little difference in the jet topology. In this scenario the 
tightest limits on SUSY particles will come from e"'"e~ colliders. 

It is well known that the non-observation of the decay of the Z° into squark-antisquark 
pairs sets a lower limit of approximately 42 GeV for the lightest squark mass . However, 
if the gluino is much lighter than Z^, a stronger bound can be put on the squark mass due 
to the effect on the Z° hadronic branching ratio of the decay 

qqg + Charge Conjugate (2) 

where the tildes indicate supersymmetric partners of the quark q and gluon g. 

The effect of squarks and gluinos on e"'"e~ annihilation into hadrons through a virtual 
photon was treated a decade ago and these results can be carried over simply to 
decay. The hadronic to electronic branching ratio of the Z^ is of the form M 



r(Z° hadrons) a3{Mz) [a3{Mz)\'^ /asiMz 
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Including the effect of virtual squark-antisquark pairs and quark-squark-gluino decay chan- 
nels, the lowest order effect of the SUSY particles is to replace as by ctefi in the linear 
term of Eq. (||) . 

asiMz) aeff(Mz) = a3(Mz)(l + d) (4) 

The 6 parameter, which is a function of the squark and gluino masses, gives the lowest 
order effect of the SUSY particles assuming the squark mass is above half the Z^ mass. 
Since the effects of SUSY particle on the other terms are not readily available we make 
the substitution of Eq. (H) in the quadratic and cubic terms of Eq. @ also, accepting 
a theoretical error of order {a3{Mz)/Tv)'^d. S can be calculated from Eqs. (3) and (8) of 
the first paper of Ref. as a function of the squark mass and the gluino mass. The 
typographical error in the last term of Eq. (8) has been corrected in our calculations. In 
these equations the quark mass is neglected by comparison to Mz and rriq. The results 
for d are presented in Fig. 1 including all quarks except the top and setting all squark 
masses equal as in Ref. H. Results for more realistic calculations in which the squark 



masses differ are discussed below. The effective strong coupling constant as deduced by 
the experimental value of the left hand side of Eq. (Q) is 

aeff(Mz) = 0.141 ±0.017 (5) 
which can be compared with the world average value 

a3(Mz) = 0.113 ±0.003 (6) 

At the one standard deviation level, Eq. (^) would imply that 

5 = 0.25 ±0.15 (7) 

One could, conservatively, take Eq. (|^ as putting an upper limit of 0.4 to the contribution 
of squarks and gluinos to S. We would not expect an appreciable fraction of this to be 
attributable to the higher order process —>■ qqgg. From Fig. 1, assuming negligible 
gluino mass, this translates into an experimental bound 

ruq > 50 GeV (8) 

Since, however, in the mi/2 — scenario for soft SUSY breaking the gluinos are 
expected to be below 1 GeV and the squarks could be between Mz/2 and Mz, it is 
interesting to consider the consequences of the assumption that squarks and gluinos are 
in fact causing the deviation of 5 from zero. Arguing against this interpretation is the 
fact that the R parameter of e"'"e~ measured at 35 GeV where there is certainly no squark 
contribution also yields an anomalously large value of as (with large errors) . We might 
suggest, however, that at 35 GeV there is room for higher twist, power law, contributions 
to the R parameter at the level of a few percent. This would be sufficient to bring the value 
at 35 GeV into line with other measurements. At the mass these should be negligible 
and therefore one might expect the hadronic decay rate of the Z^ to accurately reflect the 
perturbative QCD final states. Eq. (^, at the one standard deviation level, together with 
Fig. 1 would then imply 

50 GeV < rrig < 82 GeV (9) 

As discussed above, the uncertainties in the hadronization models built into the jet Monte 
Carlos are in the case of light gluinos such that Eq. (|^) can probably not be reliably ruled 
out by the present hadronic collider data. The sfermions, except for the stop quark, are 
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expected in a good approximation to be unmixed partners of the left and right handed 
fermions with masses 



f,RR 



f,LL 



tjiq + TTij + M|cos2/3 T3J — e/sin^6'H' , 
rriQ + rnj + M|cos2/3e/sin^6'vK- 



(10a) 
(106) 



The average squark squared mass as well as the average slepton squared mass therefore 
satisfies 



Some squark masses in the range of Eq. (P) would follow in this scenario if mo were near 
the lower end of the range in Eq. (|l|). Because of the non-linearity of 5 as a function of 
the squark masses, the effective mass in Eq. (^ for determing S lies between the lightest 
and the average squark mass. 

We ask the question: If a 5 in the range of Eq. (^ is to be produced by squarks 
and gluinos what must be the values of mo and (3 and what would be the consequences 
for slepton and top quark masses? If these consequences can be ruled out experimentally 
alternate interpretations must be sought for the 6 anomaly. Of course the anomaly is 
only a two standard deviation effect which could easily disappear with better statistics. 
In either case the parameters of the mi/2 = model would be further constrained with 
important implications for the gluino and squark masses. 

From Eq. ( pUf ) we can see that, if cos2/3 is positive, the "down-type" squarks (T3 = 
— 1/2) would be below mo and the "up-type" squarks would be above mo, with the reverse 
situation holding if cos2/3 were negative. The squared coupling to the is proportional 
to T|j -f (T3 J — 2e/sin^6'vi/)^ which is greater for "down-type" squarks than for "up-type". 
Thus, for fixed mo, a greater contribution to 6 would result if cos2/9 were positive. However, 
current experimental constraints are such that lower mo is allowed if cos2/3 is negative so 
that the largest effect consistent with current data is obtained if cos2/9 is negative and is 
dominated by the effect of "up-type" squarks. In Figs. 2 and 3 we show the d contours as 
functions of mo and the gluino mass in the cases of positive and negative cos2/5 respectively. 
However, apart from the window at very low gluino mass, most of the parameter space of 
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In the mi/2 = scenario the angle (3, the arc-tangent of the ratio of Higgs vevs is con- 
strained by current data to lie in one of two ranges [0 



-0.674 < COS2/3 < -0.385 or 0.406 < cos2/3 < 0.536 



(12) 
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Figs. 2 and 3 is ruled out by collider and beam dump searches for SUSY particles. In Fig. 4 
we show the S contours as a function of mo and cos2/3 assuming negligible gluino mass 
as suggested in the second article of Ref. . The currently allowed regions of parameter 
space in the mi/2 = model with mo < 95 GeV lie in the trapezoid in the upper right of 
Fig. 4 and in the pentagonal region in the lower right. The region to the left of the slanted 
lines of postive slope is ruled out by the current lower limit of 65 GeV on the mass of the 
selectron while the region to the left of the slanted line of negative slope is ruled out 
by the current limit on sneutrino masses. As can be seen from Fig. 4, squarks below the 
can account for the anomalously large apparent value of a^iMz) seen in decay if 
mo is in the range 

60 GeV < mo < 85 GeV (13) 

with values above 5 = 0.16 allowed only for negative values of cos2/?. 

Such a low value of mo would also have strong consequences for the top quark mass 
since minimal SUSY unification with gluinos below the Z^ leads to the effective SUSY 
mass 

Ms = 150 GeV X e-518.5(sin2eM.-0.2333)gl.85(a3-^(Afz)-0.113-i)_ ^^^-^ 

This equation with the world average values of as and sin^6'vK leads to Eq. (|l|). However, 
since we are exploring the consequences of light SUSY particles, for as{Mz) we use the 
tightly constrained value coming from the best fit to the quarkonia data allowing for light 
gluinos as{Mz) = 0.113 ± 0.001. Then requiring an Ms consistent with Eq. 13 leads 
to a tight prediction for siv?9w The experimental result for sin^6'vK is strongly dependent 



on the top quark mass [|T0 



sin^OwiMz) = 0.2327 ± 0.0004 + 0.0043 [l - 1 . (15) 

The quoted error here comes in equal parts from the experimental error and the uncertainty 
due to the Higgs mass which is assumed to lie between 50 GeV and 1 TeV. It is clear from 



Eqs. (|4|) and (|15D that lower values of mt are correlated with lower values of Ms- 



mt 



125 + 56.1 Inf —\ ±14 
V205/ 



GeV. (16) 



The observation of Ref. would increase the error here unless the GUT-scale particles 
are very nearly degenerate. Equating Ms with the average squark mass squared leads via 
Eq. (P to 

Ml = ml + 1608 GeV^ (17) 
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where, in calculating the average quark mass, we have used for mt a value near the current 
experimental lower limit since in order to bring mo into the range of Eq. (p!3|), one needs 
small values of the top quark mass. As mo increases through the range of Eq. (^3D, the top 
quark mass runs from 67 ± 14 GeV to 81 ± 14 GeV. These values are very close to the best 
fit values of mt arising from the the total unconstrained LEP data [^. However only the 
highest values here are consistent with the recent CDF limit [rT| mt > 91 GeV. Although 
the possibility of SUSY decays might loosen the CDF limit, their result taken at face value 
would disfavor a value of S greater than 0.1 coming from the mechanism discussed here 
which however is sufficient to bring aeff{Mz) into agreement with the world average value 
as{Mz) at the one standard deviation level. This would require a value of mo between 75 
and 85 GeV in agreement with the lower limit of Eq. (|l|). 

The following conclusions can be drawn from our analysis. The current experimental 
constraints on the = model of SUSY breaking are such that the discrepancy between 
the as{Mz) value from the width and the world average is consistent (at the one 
standard deviation level) with an explanation in terms of decay into quark-antisquark- 
gluino. This explanation is however not a prediction of the mi/2 = model since it 
requires that mo be near the bottom of the currently allowed range. If this explanation 
is in fact chosen by nature, many of the SUSY particles are below the Z^ and will be 
clearly seen at LEP II. The non-coUinear lepton pairs Z^ —>■ |U+/x~77 coming from the 
decay Z^ n^fl~^ + C.C should be seen even at LEP I when increased statistics are 
available. In addition the top quark should then be at the very low end of the currently 
allowed range and should be seen in future Fermilab experiments. 
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Figure Captions 

Fig. 1. Values of 5 as a function of squarlc and gluino mass. Contours of constant S are 
plotted with the value of 5 shown to the right of the contour. It is assumed that 
all squark masses are equal. 

Fig. 2. Values of 5 as a function of mo and the gluino mass for cos2/3 = 0.5. Contours of 
constant 6 are plotted with the value of S shown to the right of the contour. 

Fig. 3. Values of 5 as a function of mo and the gluino mass for cos2/3 = —0.5. Contours 
of constant 6 are plotted with the value of 6 shown to the right of the contour. 

Fig. 4. Values of 5 as a function of mo and cos2/3 for zero gluino mass. That part of 
the graph to the right of the dashed lines is the region allowed by the restriction 
that the squark masses are above Mz/2. The two enclosed areas are the regions 
allowed by Eq. ([1^) and the restrictions that the slepton mass must be above 65 
GeV and the sneutrino mass must be greater than 41 GeV. 
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